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An exact quantization rule for the bound states of the one- 
dimensional Schrodinger equation is presented and is generalized to the 
three-dimensional Schrodinger equation with a spherically symmetric po- 
tential. 
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In the development of quantum mechanics, the Bohr-Sommerfeld quantization rules 
of the old quantum theory [1] occupy a position intermediate between classical and 



quantum mechanics. The WKB approximation [2-4] is a method for the approximate 
treatment of the Schrodinger wave function with another quantization rule [1] : 

f XB kdx = (n + l/2)vr, k = j2fx[E - V(x)]/h, n = 0, 1, 2, . . . , (1) 

Jx A 

where xa and Xb are two turning points 

j^ 1 E = V(x A ) = V{x B ), EV{x), x a <x< x b , (2) 

Or 

and n is the number of nodes of the WKB wave function between two turning points. 
The half-integer number is the quantum correction to the Bohr-Sommerfeld result, 



which comes from the fact that the wave function in WKB analysis may extend into 
the classically forbidden region. The WKB method is precise for the harmonic oscil- 
lator potential, and is expected to be most useful in the nearly classical limit of large 
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quantum numbers, namely, it is good when n is large in comparison with unity. Various 
refinements have been developed to improve the accuracy of the WKB method [5-7] 
where the main modification was made on the phase loss at the turning points. Re- 
cently, Cao et al. presented a calculation method, called the analytic transfer matrix 
method (ATMM) [8,9], for the energy levels of bound states of the one-dimensional 
Schrodinger equation, where the phase contribution picked up by the scattered sub- 
waves was included in the new quantization rule. The accurate numerical calculation 
results for some examples were given there [8,9]. This is a prospective method for 
calculating the energy levels of bound states of the Schrodinger equation, but has to be 
developed and improved further. In this Letter we apply the matching condition of the 
logarithmic derivatives (which is a standard method in quantum mechanics) and the 
fundamental concept of the mathematical analysis to the problem of bound states of 
one-dimensional Schrodinger equation, and derive an exact quantization rule without 
any approximation. The quantization rule is rigorous and general for any Schrodinger 
equation with only one variable. It is a fundamental problem in quantum mechanics. 

The plan of this Letter is as follows. In section 2 we will develop the method of 
the numerical calculation for the energy levels of bound states of the one-dimensional 
Schrodinger equation by the matching condition of the logarithm derivatives of the 
wave function. In section 3 we present an exact quantization rule for one-dimensional 
Schrodinger equation. The quantization rule is proved without any approximation. Its 
validity can also be confirmed with some solvable examples. In section 4 the quantiza- 
tion rule is generalized to the three-dimensional Schrodinger equation with a spherically 
symmetric potential. The conclusion is given in section 5. 

II. MATCHING CONDITION OF LOGARITHM DERIVATIVES 

Consider the one-dimensional Schrodinger equation 

^(x) = - 2 ^[E-V(x)Mx), (3) 

where /z is the mass of the particle, and the potential V(x) is a piecewise continuous 



real function of x. The logarithm derivative <f>(x) of the wave function ip(x) is 

u \ 1 dip(x) ... 

ip(x) ax 

From the Schrodinger equation (3) we have 

-r<j)(x) = -k(x) 2 - <j)(x) 2 , E > V(x), (5) 

ax 

^-Mx) = k(x) 2 - cf)(x) 2 , E < V(x), (6) 

ax 

where 



k{x) = ^/2n[E-V{x)]/n, E>V(x), 

k(x) = ^2 l i[V{x)-E]/Ti ) E < V(x), (7) 

k(x) = k(x) = 0, E = V(x). 

It is obvious that <f>(x) decreases monotonically with respect to x when E > V(x), but 
not monotonically when E < V(x). Note that near a node of the wave function ip(x) in 
the region where E > V(x), <f>(x) decreases to — oo, jumps to +oo, and then, decreases 
again. 

Arbitrarily choosing an energy E as a parameter, we assume for definiteness that 
the potential V(x) satisfies 

V(x) = VjE, — oo < x < xj, 

V(x)E, xj < x < xa or xb < x < xp, 

V(x) = E, x = xa or x = xb, (8) 

V(x) < E, xa < x < xb, 

V(x) = VpE xp < x < oo. 

xa and x B , x A < xb, are called two turning points where V(xa) = V(x B ) = E. This 
potential can be easily generalized. 

The Schrodinger equation (3) is a linear differential equation of the second order, for 
which there are two independent solutions. In the region — oo < x < xj, one solution 
is divergent exponentially, and the other is physically admissible: 

iP(x) ~ e K ' x , <j>{ Xl ) = m = ftuM - E)/m. (9) 

Similarly, in the region xp < x < oo, the physically admissible solution is 



il){x) ~ e' KFX , <f)(x F ) = -k f = -J2/j(V f - E)/% < 0. (10) 



For the general potential, if E < V(±oo), both 4>( x i) an d — 4 > { x f) are positive . 

By making use of the fundamental concept of the mathematical analysis, we replace 
the continuous potential well with a stack of thin films each of which has a constant 
potential. We first divide the region xj < x < xa into n equal films with width d n , 
where E < V(x) and xa = %i + nd n . In the jth film, xi + jd n — d n < x < xj + jd n , 
V(x) is replaced with a constant potential Vj 



V, = V(xj + jd n - d n /2), K j = pn [V, - E]/h. (11) 

Solving the Schrodinger equation at this film, we obtain 

il>j(x) = A je KjX + Bje- KjX . (12) 

On two ends of the film, the logarithm derivatives tfij-i and tpj, which should match 
with the logarithm derivatives at the ends of the neighboring films, are 
1 dipj (x) 



<Pj-i 



tpj (x) dx 
1 dipj (x 



ti/n " 



ipj(x) dx 



x=x I +jd n 3 3 

J^. e Kj(xi+jdn-dn) _ B . e -K j (x I +jd n -d n ) 
— is . 3 J 

3 A-e K 3^ Xl+3dn ~ dn ^ + Be" K j( x i+3 d n-d„) ' 

x=x I +jd n -d n 3 '3 



(13) 



From the second formula of Eq. (13) we obtain 

A ^ Kj (x I+j d n -d n ) |^. _ ^_^ = B . e -K j {x I +0d n -d n ) { K . + ^,_^ _ 

Substituting it into the first formula of Eq. (13) we have 

Kj tanh (Kjdn) + <pj-i . . 

(fij = Kj— —. ^r- . (14) 

Kj + tfij-i tanh {Kjd n ) 
This is a recursive relation. If ifij-i is positive, tpj is positive, too. There is no zero 
both in the numerator and in the denominator of Eq. (14), so tpj is finite and non- 
vanishing. Since <f>(xi) is positive and known, one is able to calculate <fi n = 4>(xa) from 
</?o = <f>(xi) = Ki with the recursive relation (14) as j increases from 1 to n. 4>(xa) is 
positive, finite and non- vanishing (see Appendix A). The calculated precision depends 
on the number n of the films. In principle, one may obtain a precise <j>(xa) if n is large 
enough. 

Similar calculation can be made in the region ij? < x < oo. The recursive relation 

(14) becomes 

tpj - Kj tanh (Kjd n ) 

(fj-i = k>j— , / , x - (15) 

Kj — (fij tanh {Kjd n ) 



If tpj is negative, tpj-i is negative, finite and non-vanishing. Since (j)(xp) is negative 
and known, one is able to calculate y?o = <P{xb+) from (p n = 4>(xf) = —Kf with 
the recursive relation (15) as j decreases from n to 1. <P(xb+) is negative, finite and 
non-vanishing. In principle, one may obtain a precise 4>(xb+) if n is large enough. 
Physically, the fact that there is no zero of <f>(x) in the regions — oo < x < Xa arid 
x b < x < oo implies that in those two regions the wave function ip(x) is exponential 
decay. However, there may be a zero of 4>(x) in the other classically forbidden regions 
(see the end of this section). 

Now, we divide the region xa < x < xb into m equal films with width d m , where 
E > V(x) and xb = xa + md m . In the jth. film, xa + jd m — d m < x < xa + jd m , V(x) 
is replaced with a constant potential Vj 



Vj = V(x A + jd m - d m /2), kj = y/2fj, [E - Vj]/h. (16) 

Solving the Schrodinger equation at this film, we obtain 

ipj(x) = Cj sin (kjX + Sj) . (17) 

On two ends of the film, the logarithm derivatives <f>j-i and <f>j, which should match 
with the logarithm derivatives at the ends of the neighboring films, are 
1 dipj(x) 



y j-i 



ipj (x) dx 
1 dipj(x) 



ipj (x) dx 
From Eq. (18) we obtain 



= kj cot [kj(x A + jd m - d m ) + Sj] , 
= kj cot [kj(xA + jd m ) + Sj] . 



and 



x=XA+jd„ 



kj cot < Arctan J ] + kjd m \ , (19) 



/n . kA 

kjd m = —Arctan + Arctan — + qn, 

\<f>j-ij \<t>jj 

no zero of <f>(x) occurs in xa + jd m — d m < x < xa + jd ri 

1 a zero of (f)(x) occurs in xa + jd m — d m < x < xa + jd r , 



(20) 



where Arctan (3 denotes the principle value of the inverse tangent function: 

Arctan (3 = a, (3 — teaia, — n/2 < a < n/2. (21) 



Note that (f>(x) decreases monotonically with respect to x when EV(x) [see Eq. (5)]. 
If no zero of (f>(x) occurs in the film xa + jd m — d m < x < xa + jd m , 

Arctan f — — ] < Arctan ( — | < tt/2, 
\<f>j-ij \<t>jj 

Equation (20) holds with q = 0. If one zero of (f>(x) occurs in the film xa + jd m — d m < 

x < xa + jd m , 4>j-i > 0) an d 4>j < 0- Thus, 

Arctan — — I ~ tt/2, and Arctan I — I ~ —tt/2, 
\<Pj-iJ \<PjJ 

we have to add an additional it on the right-hand side of Eq. (20) such that its right- 
hand side is positive and equal to kjd m . Since the width d m of the film is very small, 
we do not consider the case where more than one zeroes of <p(x) occur in the film. 

Equation (19) is a recursive relation, with which one is able to calculate m = 
4>(xb—) from 0o — <P{ x a) as j increases from 1 to m. The calculated precision depends 
on the number m of the films. In principle, one may obtain a precise 4>(xb—) if vn is 
large enough. 

From the Sturm-Liouville theorem (see Appendix B), as E increases, 0(:r#— ) de- 
creases monotonically and 4>(xb+) increases monotonically. Choosing the parameter E 
by dichotomy such that 0(xb— ) matches with <J)(xb+), we obtain a bound state with 
the energy E. 

Cao et al. [8] derived the recursive relations similar to Eqs. (14), (15) and (19), 
and calculated some examples with more precise results than those obtained with the 
nonintegral Maslov index and the standard WKB method. Zhou et al. [9] calculated a 
problem with a one-dimensional symmetric double-well potential. Since the potential 
is symmetric, V(—x) = V(x), the solution ip(x) of the Schrodinger equation is even 
or odd in the spatial inversion, and then, its logarithm derivative 4>(x) is odd. If a 
nontrivial solution ip(x) is odd, ^(O) = and 4>(x) has to be infinity, 0(0+) = +oo. 
If a nontrivial solution ip(x) is even, -0(0) ^ and <f>(x) has to be zero because <f>(x) 
is odd, namely, 0(0+) = 0. One can calculate the energy levels of this system by the 
above method in the half space < x < oo. The different boundary conditions of 
0(0+) will give different solutions with different energies, as calculated in [9]. 



III. QUANTIZATION RULE 



In the preceding section, we divided the region xa < x < xb into to equal films, 
where E > V(x), and obtained Eq. (20) for kjd m . Summing up Eq. (20) from j = 1 
to j = to, we obtain 



7 , ^j"r, 



A^7r — Arctan [ — 

V0o, 

n'r. 



fkA A /fc 2 N 

Arctan — — Arctan — 



— h . . . — Arctan 



v m— 1 



> 



Arctan 



m— 2 , 



"m— 1 



'm— 1 



Arctan 



h 



Arctan — 

V02, 



'm— 1 



Arctan I — ^ 






(22) 



where 0o — ^(^a), 0m — 4>{xb), and A is the number of zeroes of the logarithm 

derivative <f>{x) in the region i^ < i< ig. When to goes to infinity, d m tends to zero, 

and the sum in Eq. (22) becomes an integral. Thus, we obtain a new quantization 

rule: 

f k \ 



x b 1 ( k\ 

k(x)dx = Nn + lim < —Arctan — — - 

x A ' m ^°° { \</>{ x A), 



Arctan 



\4>(%b), 



m—l 



+ E 



Arctan — 



Arctan — — 



(23) 



where <pj is calculated recursively with Eq. (19). The first term of the right-hand 
side of Eq. (23) comes from the zeroes of the logarithm derivative <f)(x) in the region 
xa < x < xb- Since 0(^)0 and </>(xb) < 0, the second and the third terms are 
vanishing as to goes to infinity if the potential is continuous at the turning points. The 
last sum denotes the phase contribution devoted by the scattered subwaves. 

The formula (23) has another expression. If one changes Sj = <5' + n/2 in Eq. (17), 
Eq. (18) becomes 



4>j-i = —kj tan kj(xA + jd m — d m ) + 5'j 
<pj = —kj tan kj(xA + jd m ) + 5'j 

Then, equations (20) and (23) become 

X 0i-i X / N 



hd 



Arctan 



kj 



Arctan — ^ + q'ir, 
\ k j) 



no zero of ip(x) occurs in xa + jd m — d m < x < xa+ jd m , 

1 a zero of ip(x) occurs in xa + jd m — d m < x < xa + jd m , 



r k{x)dx = iV'vr + lim { A rctan (^) - Arctan (^) 
+ > Arctan - r - L - - Arctan -f^- > , 

where JV' denotes the number of nodes of the wave function ip(x) in the region xa < x < 
x B . If the potential is continuous at the turning points, due to 0(^)0 and 4>{xb) < 0, 
the second and the third terms are 

lim Arctan (^A) _ Urn Arctan (^1) = W2 - (-vr/2) = vr. (25) 

m^oo \ kl J m ^°° \ k m ) 

Since 0(x) decreases monotonically in the region xa < x < xb, N = N' + 1. 

The sum on the right-hand side of Eq. (23) can be transformed into an integral 
expression: 



XB 1 I \ 1 nr I- f A ( ^1 \ . ( k 



I k(x)dx = Nir + lim < —Arctan — - — - + Arctan 

Jxa m ^°° \4>(Xa) 



111 



'*a — y \4>(xa)J \4>( x b)J j (2Q) 

[ x b <f)(x)(dk(x)/dx) y ' 

Jx A 4>{x) 2 + k{x) 2 

Two terms in the curly brackets are vanishing as m goes to infinity if the potential is 
continuous at the turning points. The sum in Eq. (24) can also be transformed into 
an integral expression. 

The quantization rule (26) is proved without any approximation, so that it is exact. 
Its validity can also be confirmed by comparing it with the following solvable exam- 
ples. Both Eqs. (23) and (24) are the formulas of numerical calculation for Eq. (26). 
Cao et al. [8,9] presented an expression similar to Eq. (24) and demonstrated it to 
be very effective in numerical calculation through two examples: the one-dimensional 
Schrodinger equation with a power-law potential [8] and a symmetric double-well po- 
tential [9]. 

Ex. 1. The harmonic oscillator potential. 

The WKB method is precise for the harmonic oscillator potential. Now, we are 
going to check our new quantization rule (26) for the harmonic oscillator potential 
V(x) = fitu 2 x 2 /2. Let 

a = \ -7T , £ = ®x, (27) 

V h 



we have 

i/> n (x) = N n e-t 2 / 2 H n (0, E n = huo{n + 1/2) , 

<p n (x) = -a£ + 2naH n _ l {£)/H n {£), 
where N n is the normalization factor, H n (£) denote the nth Hermitian polynomial. In 



the region xa < x < Xb, where — xa = Xb = \Z2rT+T/a, we have 

fc(»)(a;) = apn + 1 - £ 2 , l ' = -a A x/k {n \x). (29) 

(JjJb 

After calculation for < n < 10, we know that <p n (xA) = — re (^s)O, <f> n (x) has 
n + 1 zeroes in the region xb < x < Xb, and Eq. (26) becomes 



X A 



k {n \x)dx = (n+l)7r + $„, 
fx B (j)(x) dk^ n \x)/dx 



*" /^ 0(x)2 + A;W(a;) 2 



dx 



o [2n + 1 - e] ^n(0 2 + [2nff„_i(0 " WOf V2n + 1 - ? 

n p/R+i (2n - C 2 ) dj 
- Z J y/2n + l-C 2 n ' 

At least for < n < 10 we obtain $ n = —tt/2 by Mathematica. In fact, the first term 
in $ n is (2n — 1)tt/2, and through a variable transformation: 



2n + 1 - e = t u + v/^m , ^ 



4V2n + ltdt 
(l+£ 2 ) 2 ' 

the second integral F n is calculated to be — mi: 

F = 0, 

"^3 _2di r 1 (l + t 2 )dt 



1 Vo ,/3^F(H-£ 2 ) y 



-7T, 



Fa = 2 



F 3 = 2 



v^^F(l + e 2 ) ^o l-t 2 + £ 4 

v^ -4(4£ 2 + 5)d£ 
o v/5^F(4£ 4 + 4£ 2 + 5) 

_ yi (l+t 2 )(5-6t 2 + 5t 4 )cit = _^ 
7o 25-76£ 2 + 118£ 4 -76£ 6 + 25£ 8 
v^ -6(4£ 4 + 6£ 2 + 9)d£ 

y 1 (1 + £ 2 )(247 - 748£ 2 + 1146t 4 - 748t 6 + 247£ 8 )d£ 

Jo 1444 - 8052t 2 + 20652t 4 - 27512t 6 + 20652t 8 - 8052t 10 + 1444t 12 

-3tt. 



Thus, we demonstrate that the quantization rule (26) is the same as Eq. (1) for the 
harmonic oscillator potential. 

Ex. 2. The square well potential. 

Discuss a finite square well potential V(x) 



V(x) 



Va x < —n, 

V B X>7X, (31) 

—71 < X < IX. 



The logarithm derivatives at the turning points xa = — 7r and xb = tt are 



4>{x A ) = K T = pn(V A - E)/h, 4>{x B ) = -k f = -pn(V B - E)/h, (32) 

when E < Va and E < Vb- The solution to the Schrodinger equation is 

if> n (x) = sin (k {n) x + 5 {n) ) , (33) 

and its logarithm derivative <f> n (x) is 

<f> n (x) = k {n) cot (k {n) x + 5 {n) ) . (34) 

k^ n > and 5^ n > are determined by the matching conditions at the turning points x a and 

x B , 

tan (-k^ix + 5 {n) ) = k {n) /K U tan (h {n) TX + 5 {n) ) = -k {n) /K F . 

Hence, we obtain 

£>) = - {Arctan (k {n) / k/) + Arctan (k {n) / k f ) } /(2tt) + n/2. (35) 

4> n (x) has n zeroes in the region — ix < x < tx, and k(x) takes the constant value k^ 
in the region, dk(x)/dx = 0. The energy levels are E n = (hk^j /(2/x). In terms of 
Eq. (35), the right-hand side of Eq. (26) is calculated to be 



nix — Arctan — — - + Arctan - 



k 



M x a)J \(f>(x B )J (36) 

= mx - {Arctan (jfeW/wj) + Arctan (Jc^/kf)} = 2ixk^ . 

Due to the constant k(x), the left-hand side of Eq. (26) is equal to the same value: 

k{x)dx = 2ixk (n) . (37) 



X B 

XA 



This quantization rule is different from both the Bohr-Sommerfeld one and that given 
by the WKB approximation. When Va and Vb tend to infinity, due to Eqs. (32) 
and (35) k^ goes to n/2, and the quantization rule (37) coincides with the Bohr- 
Sommerfeld one. 

IV. THREE-DIMENSIONAL SCHRODINGER EQUATION 



Consider the three-dimensional Schrodinger equation with a spherically symmetric 
potential. After separation of the angular part of the wave function, 



^r)=r- 1 R(r)Yi(9, V ), 

the radial equation of the Schrodinger equation is 
d 2 R{r) 



dr 2 



| {£ - U(r)} R(r), U(r) = ^ + 1} + V(r). 



(3* 



(39) 



Since Eq. (39) is similar to Eq. (3), its energy levels can be calculated by the matching 
conditions of the logarithm derivatives, where the logarithm derivative is defined as 

_ x rfi?(r) 



(r) = R(rY 



dr 



(40) 



As an example, we discuss the problem of the hydrogen atom, where the potential 

V(r) is 

V(r) = --. (41) 

r 

When r — > oo, we have 



R(r) ~ exp{-r^\E\/h}, 0(r) ~ -^\E\/h. 



When r — > 0, we have 

R(r) ~ r t+1 



fxe 



h z u+i) 



(r)~(£+l)/r. 



By the method of matching condition of logarithm derivatives, one is able to calculate 
the energy E with Eqs. (14), (15) and (19). 

On the other hand, if the solution to Eq. (39) has known, we are able to check 
whether the quantization rule (26) holds. For the energy E n 



E n 



/j,e 

2h 2 n 2 ' 



n = 1,2,3,..., 



(42) 



we have the solution [1] 



Rni{r) = N ne e-^ 2 p E+1 L 2 n %\p), p 



2pe 2 



2l+\ 



(43) 



where N n e is the normalization factor and L n + e (p) is the associated Laguerre polyno- 
mials. 

When £0, the turning points r A and r B satisfying U(r A ) = U(r B ) = E n are 



ta 



r B 



nh 2 

nh 2 
2pe 



PA, 



;Pb, 



p A = 2\n- 



n 2 -£(£+!) 



p B = 2{n+ n 2 -£{£+!) 



1/2 
1/2 



(44) 



When £ = 0, we define r A = with U(r A ) = — e 2 /r A ~ — oo. r B with U(r B ) = E n still 
satisfies Eq. (37). The momentum k n i(r) is 

(45) 



J 1 Icy 

k ne (r) = 27 {(P - Pa) (pb - p)} ' 



In the Schrodinger equation for the hydrogen atom, the quantization rule (26) becomes 

' ' M + tat-f^)! 

(46) 



k n Ar)dr = Nn + lim < — Arctan . 
>:\ ' m - co { \Mr A ) / 

TB Mr) (dk ne (r)/dr) 



s (p n e(r B ) 



-dr. 



'r A <j>ne(r) 2 + k nl [rf 
Calculating the left-hand side of Eq. (51), we obtain with Eq. (50) 

r-r B 



k n e(r)dr = — {(p - p A ) (p B - p)} 1/2 

r A J PA ^P 



n 



v. 



7T. 



(47) 



Now, we calculate the right-hand side of Eq. (51). The number iV of zeroes of 4> n i(r) 
in the region r A < r < r B is (n — £). The logarithm derivatives at the turning points r A 
and r B are non-vanishing. When £ = 0, r A is not the turning point, but <f) n o(r A ) = 1/r. 
Thus, we have 



lim Arctan 



lim Arctan 



lim Arctan 




0, 
0, 



lim Arctan <j r\J2p [E n + e 2 /r)/Ti \ = 0, 



>0, 



0. 



(4£ 



s (j>no(r A ) i 

Calculating the integral in Eq. (51) by Mathematica, at least for n — 1, 2, 3, and 
£ < n, we obtain: 

rrB <Pne{r) {dk ne (r)/dr) 



-dr 



£-^£{£ + T) 



'ta <j) n e(r) 2 + k ne (r) 2 
Therefore, the quantization rule (46) holds for the hydrogen atom 



7T. 



(49) 



V. CONCLUSIONS 

In this Letter, with the matching condition of the logarithm derivatives and the 
fundamental concept of the mathematical analysis, we proposed a formula (23) for nu- 
merically calculating the energy levels of bound states of the Schrodinger equation in 
one dimension. Calculating the integral form of Eq. (23), we obtained an exact quan- 
tization rule (26) for the bound states of the one-dimensional Schrodinger equation. 
The exact quantization rule was generalized to bound states of the three-dimensional 
Schrodinger equation with a spherically symmetric potential. The quantization rule was 
confirmed by checking some examples where the solutions of the Schrodinger equation 
are known. Two examples of numerical calculation for the one-dimensional Schrodinger 
equation with a power-law potential [8] and with a symmetric double-well potential [9] 
demonstrated that the exact quantization rule is very effective in numerical calculation. 

Cao et al. [8,9] derived a formula similar to Eq. (24). However, their formulation 
contained some unclear points. The integer N in Eq. (29) of Ref. [8] is unclarified. 
They introduced the "exponentially decaying coefficients" Pj, but the physical meaning 
of Pj is rather ambiguous. In fact, Pj are nothing but the logarithm derivatives of the 
wave function (in their region of < x < xc) or that multiplied with —1 (in their 
region of xd < x < xs)- This is the reason why the coefficients should be matched at 
the turning points, which was not explained clearly in Ref. [8]. Finally we would like 
to point out that the series form (23) is an approximate formula of the integral form 
(26) derived in the present Letter, which is the exact quantization rule. 
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Appendix A Property of <J)(xa) 

From the recursive relation (14) and 0(x/)O, we have proved that 4>(xa) is non- 
negative and finite. Now, we are going to prove that 4>(xa) ^ by reduction to 
absurdity. When x < xa and x near xa, consider the leading term in the power series 
in (xa — x). Due to k(xa) = 0, we have k(x) ~ C(xa — x) a , where CO and aO. If 
4>{xa) = 0, from Eq. (6) we have <J){xa) — —C 2 (xa — x) 2a+1 /(2a + 1) < 0. It conflicts 
to the fact that 4>(xa) is non-negative. 

Appendix B The Sturm-Liouville Theorem 

Denote by ip(E,x) the solution of the Schrodinger equation (3) with the energy E. 
Multiplying Eq. (3) with ip(Ei,x) we have 

d 2 2u 

^(E u x)^(E, x) = ~^[E- V(x)] if,(E, x)^(E 1 , x), (Bl) 

Exchanging E\ and E and subtracting from it by Eq. (Bl), we obtain 

(B2) 
When E < V(— oo), the boundary condition gives that both solutions ip(E,x) and 
ip(Ei,x) are vanishing at negative infinity. Integrating Eq. (B2) from — oo to x^ — , we 

obtain 

1 \ , d^(E u x) ME,x)\ 

'ip(E,x) w(E u x) 



E] — E dx dx 

2/x [ x b 



n 



, j ip(E,x)ip(Ei,x)dx. 



Taking the limit as E\ goes to E, we have 



d<f>(E,x) 



dE 



» J. 



di[)(E,x) 



dE\ip(E,x) dx 
2/j, 



hy(E : x B ) 2 J~oo 



[ XB ip{E,x) 2 dx < 0. 

J — oo 



x=a;s- 



(53) 



Namely, at a given point x B — , the logarithm derivative (f)(E,x) of the wave function 
ij)(E,x) decreases monotonically as E increases. 

Similarly, when E < V(oo), the boundary condition gives that both solutions 
ijj(E,x) and if)(Ei,x) are vanishing at positive infinity. Integrating Eq. (B2) from 
Xb+ to oo, we obtain 



1 



dip(E!,x) dip(E,x) 

ip{E, x) ip{E 1: x) 



Ei — E [ T dx dx 

Taking the limit as E\ goes to E, we have 



X=X B + 



—j / tp(E,x)tp(E 1 ,x)dx. 



d(f)(E,x) 



dE 



2^ 



ip(E,x) 2 dxO. 



{BA) 



X=XB+ htlj(E,x B ) 2 Jx B 

Namely, at a given point x B +, the logarithm derivative <j)(E,x) of the wave function 
i/j(E,x) increases monotonically as E increases. 



